Excess quantum noise due to mode nonorthogonality in dielectric microresonators 
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This work presents a theory of the frequency-resolved light emission of active two-dimensional 
dielectric microresonators, which are characterized by a highly non-paraxial mode structure and 
frequently feature a position-dependent dielectric constant and non-uniform gain. The Lorentzian 
intensity profile is characterized by an appropriately generalized Petermann factor, a renormalized 
peak position, and the cold-cavity resonance lifetime Y. The theory also delivers a relation of Y 
to the laser threshold that improves earlier phenomenological expressions even for the case of a 
homogeneous medium. 
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Breit-Wigner theory is a cornerstone of the study of 
open quantum systems and finds numerous applications 
from atomic and nuclear physics to mesoscopic transport 
and optical or acoustical resonators. This perturbative 
theory starts from the isolated situation and delivers a re- 
lation between the total scattering cross section and the 
resonance linewidth Auj. Because of its perturbative na- 
ture, Breit-Wigner theory breaks down for moderately 
open systems. A context which is particularly suited 
to explore the ensuing new physics are lasers, which are 
known to emit light around a sharply defined frequency 
ou with a Lorentzian intensity profile, 
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Here I = J duoI{uj) is the total output intensity. The 
linewidth is ultimately limited by the spontaneous emis- 
sion of photons. For almost lossless resonators, the 
Schawlow-Townes formula relates the quantum-limited 
linewidth Ac^st = r 2 /27 to the total output intensity / 
and the cold-cavity decay rate Y [1]. This prediction is 
based on Breit-Wigner theory, and amounts to one spon- 
taneously emitted noise photon per cavity mode. The 
explicit factor of 1/2 accounts for the amplitude suppres- 
sion of field fluctuations due to the active feedback of the 
medium, and is absent below the lasing threshold. 

That the sensitivity to quantum noise increases for 
more open resonators was first realized by Petermann 
@], who considered gain-guided semiconductor lasers. He 
predicted that the linewidth is enhanced by a factor of 
K > 1 such that Auj = KY 2 /2I. Siegman [3] developed 
a general framework to study open-sided resonators with 
a paraxial mode structure and arrived at the expression 
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which relates the Petermann factor to the inverse con- 
dition number of the transverse resonance wavefunction 
tpo. This establishes a direct connection to a mathemat- 
ical measure of mode nonorthogonality. Subsequently, 
several groups demonstrated the excess noise in experi- 
ments on various open-sided resonator geometries with 
different cross-section 0, H, @, 0, E[ • 
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FIG. 1: (Color) Inset: the white circles indicate the dielectric 
interfaces of an annular resonator with radii R2 — 0.6 Ri and 
eccentricity d = 0.22 Ri . The refractive index in the inner disk 
is 712 = 3.3, while in the annular region Rem is fixed to 1.8. 
The color-coded density plot shows a TM-polarized resonator 
mode at its threshold (m,o = 1.8 — i0.029, ujqRi — 6.45). 
Main panel: frequency-resolved intensity close to the thresh- 
old (suitably scaled to be independent of the small detuning 
of Imni). The exact numerical result (solid black curve) is 
a Lorentzian whose width and height can be accurately de- 
scribed using the theory developed in this paper (solid red 
curve on top of the black curve) [Eq. J3} with parameters de- 
termined from Eqs. J4|-(j7| (TM polarization)]. A discrepancy 
in the height is observed when the Petermann factor is deter- 
mined from the homogeneous expression K , Eq. © (dashed 
blue curved), or when it is completely neglected (dotted blue 
curve, corresponding to the Schawlow-Townes result). 



The past decade has seen the development of a new 
class of microlasers, in which the confinement is due 
to internal reflection at dielectric interfaces [1, El, EH 
El, El, EI El, EH- The typical feature size s is of- 
ten not much larger than the optical wavelength A. 
Due to their shape, these resonators support effectively 
two-dimensional, highly non-uniform, anisotropic modes 
which are of a decisively non-paraxial character. These 
systems are currently under intense experimental and 
theoretical investigations. Until very recently, how- 
ever, the consequences of mode nonorthogonality have 



2 



bee n ig nored. Models based on random-matrix theory 
[l7l Il8l [l9|, [2(j only partially account for the structured 
wave patterns of typical microresonators. The first steps 
into this direction have been taken in recent numerical 
investigations of stadium billiards [21] and spirals 12211 . 
which focus on the case of two cross-talking modes [23| 
and demonstrate that the Petermann factor can be arbi- 
trarily large even for high-Q modes. Experiments on this 
subject have just picked up [24|. Therefore, there is a de- 
mand for a systematic theory of mode nonorthogonality 
and its consequences for the quantum noise that accounts 
for the specific features of dielectric microresonators. 

In this paper I fill this gap and develop a theory of the 
frequency-resolved light emission that applies to general, 
nonparaxial modes in such resonators and furthermore 
admits spatial non-uniformity of the dielectric constant 
and the gain (I assume single mode lasing and there- 
fore also ignore degeneracies of lasing modes). The re- 
sult entails that formally, the concept of the Petermann 
factor emerges only in the limit of a small wavelength 
(or a large resonator), A/s C 1. In reality, however, the 
relations given below already work for relatively small 
resonators (for illustration see Fig. [TJ which presents re- 
sults for an annular microresonator where both the di- 
electric constant as well as the gain are nonuniform.) The 
Lorentzian intensity profile can then be parameterized by 
a few characteristic parameters, 
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This result contains an appropriately defined Petermann 
factor in the form of a generalized condition number, 
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(TM polarization), (4a) 



(TE polarization), (4b) 



where no is the complex refractive index at threshold, 
and ^0 is the associated resonance wavefunction. The 
resonance frequency at threshold is denoted by c^o- Below 
threshold, the linewidth is given by 
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where n is the refractive index at the working point. 
(Throughout this work the vacuum speed of light c = 1; 
above threshold, the linewidth is further reduced by a 
factor of 1/2 [25], as in the Schawlow-Townes case.) The 
theory also predicts a systematic line shift 
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which is of the same order as the linewidth and generally 
only disappears for large resonators filled with a homoge- 
neous medium. The cold-cavity characteristics T and Suo c 
are obtained from the relation Soo c — iT/2 — 5Q C where 
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As a notable side product, the latter expressions general- 
ize and improve the common phenomenological relation 
Imno = — (2T/uo) Re no (TM) between the cold-cavity 
lifetime and the laser threshold of homogenous systems. 

All considerations in this paper are based on the 
effective-medium approach, where the geometric and 
material-specific properties of a dielectric microresonator 
are described by a position-dependent refractive index 
n(r). Let us assume n(r) = 1 outside the resonator. In- 
side the resonator the refractive index can be complex, 
with Imn < in the amplifying (active) regions of the 
medium. For an effectively two-dimensional resonator, 
the classical electromagnetic field is represented by a 
scalar wavefunction ip(r]uo) which fulfills the Helmholtz 
equation [26[ 
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(TM), (8a) 
C = Vn- 2 (r)V + cj 2 (TE). (8b) 

Consider the case that Ren(r) is fixed while Imn(r) 
can be controlled via pumping. At the laser threshold 
Imn(r) = Imno(r), Eq. (|HJ) admits a solution ^o satis- 
fying purely outgoing boundary conditions with real an- 
gular frequency uoq. Below the laser threshold, the res- 
onance frequency Vt associated to this solution becomes 
complex, where — 2Im^ = 7 is the decay rate. The 
cold-cavity decay rate Y is obtained for Imn = (the 
complex cold-cavity resonance frequency is denoted by 
tt c =u c -iT/2). 

Quantum optics adds spontaneously emitted photons 
to this picture, which are generated and amplified by 
the active medium. Even in absence of external illumi- 
nation, the resonator then emits photons of frequency 
co w Re £7. A convenient quantum-optical framework to 
study this radiation is provided by the input-output for- 
malism, which delivers the output intensity [TH, [33] 
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in terms of the scattering matrix S(uo) of the classical field 
(assuming complete population inversion in the medium). 

For the cold cavity the scattering matrix is unitary, and 
the passive resonator does not emit any radiation. For 
a finite gain the scattering matrix departs from unitar- 
ity because the particle flux is no longer conserved due to 
the stimulated emission of photons in the active medium. 
Even if the cold-cavity resonances are strongly overlap- 
ping for the cold cavity (as is typical for low-Q modes 
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with T larger than the mode spacing), well-resolved res- 
onances appear when one steers the system close to the 
laser threshold. Dielectric microresonator design aims 
to support high-Q modes, which can be resolved even 
deep below threshold [U, [13]; in that case, however, the 
threshold is reached more quickly than in the case of low- 
Q modes. Both cases can therefore be approached by 
a systematic expansion around the threshold condition 
[28[, which delivers the output intensity 
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This expression formally diverges at the complex fre- 
quency Q = coq + SQ where 
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The output intensity is hence a Lorentzian [Eq. JTJ)] 
whose center uj = Re £7 is displaced from uoq by a sys- 
tematic line shift 5uo = ReJ^l [see Eq. (|6j)], which is 
generally of the same order as the width Auj = — 2 ImJfi 
of the Lorentzian [see Eq. J5])]. The latter can be related 
to the total intensity I = J I(uo)duo via the relation 
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This expression contains integrals similar to those in the 
expression for Kq [Eq. JS|)], but weighted by the refractive 
index. Moreover, instead of the transverse mode profile, 
now represents the resonance wavefunction in the res- 
onator plane, which encodes the full, non-paraxial mode 
structure. At this point, the integral in the denominator 
is not restricted to the interior of the resonator, and the 
result does not feature the cold-cavity decay rate. For 
the latter reason Eq. ([T2]) cannot be used to extract a 
generalized Petermann factor. 

In order to make contact with the conventional theory 
of the linewidth we now consider the regime A <C s of a 
large resonator. Under this condition, typical resonances 
have a small cold-cavity decay rate, Y <C c^o, and the laser 
threshold is attained at a small gain, |Imno| <^ Reno- 
Equation ([TT]) then applies all the way down to the cold- 
cavity limit, and can be further recast to deliver the 
complex cold-cavity resonance frequency Vt c = 00$ + SQ C , 
where SQ C is given by Eq. ([7j). This result is perturba- 
tive in the gain but nonperturbative in the openness of 



the system, which is encoded in cj ? and i/jq. Note 
that in Eq. Q, the numerator is determined by the ac- 
tive region, while the denominator also depends on the 
passive regions — including the region outside the res- 
onator, where integrals are well defined because n is 
real [29|. For homogeneous systems, the simple relation 
Q c = uo + i( Imno / Re no)cc;o (TM) is often applied (this 
also entails the relation for the laser threshold discussed 
above). In comparison to Eq. ([7j), this simple expression 
is only valid when one can ignore to the exterior contri- 
bution to the denominator, which in practice requires a 
very large system size. 

With Eq. 0, linewidth ([12]) simplifies to 
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which features the cold-cavity decay rate Y = — 2Imc)^ c 
along with the cold-cavity line shift 5cj c = Ke5Q c . (Far 
above the laser threshold, the linewidth is again halved 
because amplitude fluctuations are suppressed by the 
nonlinear feedback with the medium [25]). The Peter- 
mann factor, hence, takes the form of Eq. For TM 
polarization, the conventional form [Eq. (|2j)] of the Pe- 
termann factor is recovered when the active medium fills 
the resonator homogeneously (so that no does not de- 
pend on position). This conceptual relation hence applies 
even when the mode structure in the resonator plane is 
strongly non-paraxial. This is also the case for TE po- 
larization or an inhomogeneous medium, where Eq. (j4]) 
relates the excess noise to an appropriately generalized 
measure of mode nonorthogonality (satisfying K — > 1 in 
the limit of a closed resonator). For an inhomogeneous 
medium, Eq. (j4]) shows that the excess noise is generated 
in the active regions (where Imno 7^ 0). 

Model system. In order to assess the predictive power 
of the above analysis I turn to a paradigmatic example 
of a microresonator with a nonuniform refractive index, 
which at the same time displays a rich variety of res- 
onator modes. This is provided by the annular microres- 
onator geometry, defined by two eccentric circular inter- 
faces such as indicated in the inset of Fig. [TJ In order 
to exploit the full generality of the formalism presented 
above, I assume that the gain is constrained to the annu- 
lar region between these interfaces (Imni = —in), but 
vanishes in the interior disk (Imn2 = 0), and further- 
more set Reni = 1.8 and Ren2 = 3.3. The radii of the 
two circles are related by R2 = 0.6 Ri, and the eccen- 
tricity is d = 0.22 R\. The finite eccentricity breaks the 
rotational symmetry and hence leads to nonintegrable 
classical ray dynamics. The classical phase space ac- 
commodates domains of stability embedded into regions 
of chaotic instability and, hence, exhibits the full com- 
plexity of generic dynamical systems. The corresponding 
wave dynamics reflect this complexity in a rich set of res- 
onance wavefunctions, which also includes examples with 
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a highly directional far-field radiation pattern [33, HH . 

Numerical techniques make it possible to investigate 
a large number of modes, especially in the interesting 
regime where o;o^i is moderately large [30]. Here I focus 
on a representative example of a resonator mode with TM 
polarization, which reaches its threshold at no = 0.029, 
cooRi = 6.45. In the main panel of Fig.[TJ the solid black 
curve shows the intensity profile close to threshold (for a 
small detuning of n) , which is directly obtained from Eq. 
([9]) via a numerical computation of the scattering matrix. 
The curve is scaled so that it becomes independent of the 
detuning when n — n is small. The red curve is obtained 
from the theory developed in this work: a Lorentzian of 
form (|3]), where the parameters are determined by Eqs. 
©-([7]) (for TM polarization). This curve lies on top of 
the black curve so that the latter is barely visible. The 
theory of this paper hence applies even though the res- 
onator is only moderately larger than the wavelength. 
The other curves illustrate the necessity to adopt the re- 
sults developed in this paper. The dashed blue curve is 
obtained when the Petermann factor is taken of the ho- 
mogeneous form Ko [Eq. ([2])], while the dotted blue curve 
shows the result when mode nonorthogonality is entirely 
ignored (as in Breit-Wigner theory, which delivers the 
Shawlow-Townes formula) [32| . 

For completeness I remark that in the case of a homo- 



geneous circular disk resonator (corresponding to R2 = 
0), Eq. ([2]) applies. In polar coordinates, the radial de- 
pendence is given by Bessel functions, while the con- 
straint on time-reversal symmetry singles out modes with 
a standing- wave angular dependence sin(ra(0 — 0o)) [30[. 
For TM polarization and uoR ^> 1, Petermann factor ([2]) 
then takes the asymptotic form K = smh 2 x/x 2 where 
x = 21m (nuoRi) and therefore approaches unity for 
modes with a large Q factor. 

In summary, I have explored the consequences of mode 
nonorthogonality in dielectric microresonators for the 
general case of a highly non-paraxial mode structure 
and nonuniform material properties. This theory leads 
to an appropriately generalized Petermann factor [Eq. 
(|2|)]. Compact expressions ©-(0) can be used to inform 
the interpretation of experiments beyond the simple ef- 
fective phenomenological models employed to date. An 
open question (beyond the illustrative example of Fig. [I]) 
that can be pursued on basis of these expressions is how 
complex quantum dynamics (regular versus chaotic wave 
patterns) express themselves in the mode nonorthogonal- 
ity. 

This work was supported by the European Commission 
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APPENDIX A: FURTHER DETAILS OF THE 
DERIVATION 

This appendix details some technical steps in the 
derivation of Eq. ([TP)) . The analysis exploits that near 
resonance, the outgoing radiation is dominated by the 
resonance wavefunctions, i/j « ^o- To be specific, let us 
consider a disk C of radius R — > oo that completely covers 
the resonator, so that outside the disk n(r) = 1. Using 
polar coordinates r, tp in a coordinate system with origin 
at the center of the disk, a basis of scattering wavefunc- 
tions is then given by 



(Al) 



where H is the Hankel function and m = 0, ±1, ±2, . . .. 
The index a takes the value a = + for outgoing scattering 
states while a = — for incoming scattering states. 

The expansion of the scattering wavefunction in the 
basis (jAip defines, respectively, the in- and outgoing scat- 
tering amplitudes a m and b m: 



^(r) = ^[6 m 0W(r,(^)+a m ^(r,^) 



(A2) 



Adopting a vector notation the amplitudes are related 
via the scattering matrix 5, b = 5a, which is found by 
matching the scattering states to solutions of Eq. (|8a|) in- 
side the disk. The expansion coefficients of the resonance 
wavefunction in the scattering region are denoted by bo. 

The following considerations rest on a consequence of 
the divergence theorem, 
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which holds for arbitrary well-behaved functions ^2, 
and the orthogonality property 
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of the scattering basis functions (<9C denotes the circu- 
lar boundary of the disk C, which eventually moves to 
infinity) . 

Due to time-reversal symmetry, the scattering matrix 
is symmetric, S = S T . To see that this symmetry holds 
for complex n, consider two arbitrary solutions ipi, ip2 
of the Helmholtz equation, and apply the definition of 
the scattering matrix as well as Eqs. (|A3p and (|A4ft to 
evaluate 

- 8za^(S - S T )ai = ^ (^2^1 - ^lAta) dr = 0. (A5) 

Close the resonance, the outgoing radiation will be 
dominated by the characteristics of the resonant state. 
Hence, S ~ abobj, where a diverges at the resonance so 
that other outgoing wave components can be ignored. 



The evaluation of the coefficient a requires to compute 
the overlap between the resonance wavefunction ipo and 
an arbitrary near-resonant solution ip « ^0, which also 
possesses a small incoming component. Assume that i/j 
is normalized so that the dominant outgoing radiation is 
described by the amplitudes bo. The scattering matrix 
then takes the form S ~ bob r /(b r a), while the intensity 
becomes 
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In both expressions the denominator is independent of 
the choice of the near-resonant wavefunction ip: 



[M^)-^(^o)]dr 

ip 2 (n 2 u 2 - nlul ) dr, ( A7) 



where in the last step we used the Helmholtz equa- 
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tion. As we are close to resonance, 
2tiqUJ 2 (uj/ujo + n/no — 2). 

We now once more exploit time-reversal symmetry, this 
time in the form that i/jq provides a purely incoming so- 
lution of the Helmholtz equation where no 
fore, 



0^0 



There- 
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[^o(A^)-^(A^ )]dr 



= / \M\nl-nl 2 )uldv. (A8) 



Equation (jlQaj) now follows from Eqs. (|A6|) . (|A7|) . 
and (|A8|) . The techniques described above can also be 
adapted to modes with TE polarization which are gov- 
erned by the Helmholtz equation ([8b]) : this then results 
in Eq. (pip. 



APPENDIX B: COMPARISON TO OTHER RESONATOR MODES 
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FIG. 2: (Color) Same as Fig. \T\ but for other resonator modes of the annular resonator. All modes have even parity with 
respect to the horizontal resonator axis. The mode with uoRi — 6.45 is almost degenerate with the odd-parity mode of Fig. [1] 



